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Abstract 

We point out a precise connection between Brownian motion, Chern-Simons theory on S"^, 
and 2d Yang-Mills theory on the cylinder. The probability of reunion for N vicious walkers on 
a line gives the partition function of Chern-Simons theory on S'^ with gauge group U{N). The 
probability of starting with an equal-spacing condition and ending up with a generic configura- 
tion of movers gives the expectation value of the unknot. The probability with arbitrary initial 
and final states corresponds to the expectation value of the Hopf link. We find that the matrix 
model calculation of the partition function is nothing but the additivity law of probabilities. 
We establish a correspondence between quantities in Brownian motion and the modular S- and 
T-matrices of the WZW model at finite k and A^. Brownian motion probabilitites in the affine 
chamber of a Lie group are shown to be related to the partition function of 2d Yang-Mills on the 
cylinder. Finally, the random-turns model of discrete random walks is related to Wilson's pla- 
quette model of 2d QCD, and the latter provides an exact two-dimensional analog of the melting 
crystal corner. Brownian motion provides a useful unifying framework for understanding various 
low-dimensional gauge theories. 



1 Introduction 

Low-dimensional gauge theories are a useful playground for understanding quantum field theory. 
They also play an important role in string theory, and in particular in topological string theories. 

Relations between statistical mechanical systems and quantum field theory are well-known 
and have proved useful on both sides (see, for example, [1]). A well-known example is provided 
by conformal field theories in two dimensions. For the case of higher-dimensional topological 
theories, however, such examples have until recently remained somewhat limited. An important 
recent breakthrough was the realization [2] that the closed topological A-model vertex discovered 
in [3] is related to a melting crystal corner [4]. In fact, D-branes can be included and have a 
natural interpretation as defects [5]. In [4] this melting crystal picture was used as a definition 
of 'quantum Kahler gravity', the quantum gravitational theory describing fluctuations of the 
Kahler structure and topology while keeping the complex structure fixed. In [6], certain types of 
random walks were used to describe the combinatorics of triangulations of 2+1 quantum gravity. 
A review of the manifold connections between random walks/Brownian motion and conformal 
field theory and two-dimensional quantum gravity can be found in [7]. 

In this letter we provide further examples of precise reformulations of statistical mechanical 
systems in terms of gauge theories that are topological or close to topological. Full details will be 
given elsewhere [8]. Since the gauge theories in question have string theory realizations [9, 10], 
it is our hope that this relation can be useful in string theory. 

The first example relates Brownian motion of N vicious — that is, non- intersecting — 
walkers on a line [11] with certain boundary conditions - or, equivalently, one particle moving 
in the Weyl chamber of a simply-laced, compact Lie group G - to Chern-Simons theory [12] 
on S"^ for the corresponding group. The correspondence works for the partition function, the 
expectation value of the unknot, and the expectation value of the Hopf link invariant, the basic 
reason being the correspondence between the WZW modular S- and T- matrices and Brownian 
motion quantities. The additivity law of probabilities can be reinterpreted as a matrix model 
derivation of the partition function of Chern-Simons theory on S'^, where the repulsive force of 
the Vandermonde interaction translates into an effective repulsion exerted by the walls of the 
Weyl chamber. The correspondence is at finite values of k and N, where N is the number of 
movers and ^f^ = jrr77 ~ ~t ^^ ^^^ inverse time. This is reminiscent of relation found in [2], 
where t played the role of the temperature. 

The second example concerns a walker moving in the fundamental Weyl chamber of an affine 
Lie algebra. This can be reinterpreted as the partition function of 2d Yang- Mills on the cylinder, 
where the initial and final positions correspond to states at the two ends of the cylinder, and now 
time is related to the area. Our derivation reproduces formulas in the mathematical literature 
which make the structure of the partition function much more transparent than the usual sum 
over representations. In particular, the partition function is shown to be an affine character, 
and so the modular properties of 2d YM on the cylinder are most explicit in this formulation. 

The third example is the case of discrete random walkers rather than continuous Brownian 
motion; we analyze the random turns model [11, 13]. Its relation to randomly growing Young 
tableaux and (by Poissonization) to lattice QCD2 are known in the mathematical literature. We 
reinterpret the lattice QCD2 as a two-dimensional melting crystal, where the energy cost for 
removing a particle is the (logarithm of) the gauge coupling. 



2 Brownian motion and Chern-Simons theory 

In this paper we will be concerned with Brownian motion of N movers on a line. We will study 
so-called vicious walkers [11]. These are random walkers whose trajectories are not allowed to 
intersect. As is well-known, the intersection properties of such a random walk process play an 
important role in quantum field theory (as reviewed in detail in [7, 1]). 

Let us first consider the case of a single free walker performing Brownian motion on a line. 
The probability for it to travel from y to x in time t is given by 

Here, D is the diffusion coefficient of the medium, which from now on we will set equal to 2 • 
This probability, and its generalizations, is the basic quantity that we will study in this paper. 
It can be obtained by a variety of methods; for an overview see e. g. [1, 14]. Let us however 
mention some of its important properties. First of all, it can be obtained as the continuum 
limit of a discrete random walk, where at each tick of the clock the particle can travel a fixed 
finite distance left or right with equal probability (see also section 3). The probability is then a 
binomial distribution whose continuum limit is Gaussian. Further, pt{x,y) satisfies the diffusion 
or heat equation 

d 

— pt{x,y)=DApt{x,y), (2) 

where A is the Laplacian in x. Finally, pt{x, y) turns into a Dirac delta function of the position 
X — y when t tends to zero. 

The above has an obvious higher-dimensional generalization: 

PtA^,y) = ^^^^^^^^ e-^"^, (3) 

and N is the dimension. Equivalently, this can be regarded as the product of the probabilities 
of N single movers on a line, i.e. as the probability for A'^ movers on a line to start at positions 
yi, . . . , yjv and end up at xi, . . . , xjy after time t. 

The particular process that we will relate to Chern-Simons theory is that of N vicious walkers 
on a line. Walkers are vicious [11] if they annihilate each other when they meet. Thus, we will 
impose a non-intersecting condition and compute a probability for these walkers to walk from 
one configuration to another without ever intersecting. If we denote their coordinates by Aj, 
i = 1,...,N, they satisfy Ai > A2 > . . . > A^y. Alternatively, this process can be regarded 
as motion of a single particle in the fundamental Weyl chamber of U{N). The particle starts 
moving at position /ij satisfying /xi > ^2 > • • • > I^N, and is required to stay within the Weyl 
chamber. The process stops when the particle hits one of the walls. One then computes the 
probability of going from an initial position Hi to a final position Aj staying always within the 
chamber. This is given by [11]: 

Pt,N{X, ^) = (^2TrtW/2 e~^^^ det |e^''^^/*|i<i<j<^ . (4) 

Obviously, this probability is symmetric under interchange of initial and final boundary condi- 
tions. 

Let us now evaluate this amplitude in a specific case. We take the same initial and final 
boundary conditions, i.e. fi = X, and further an equal spacing condition, that is, Aqj — Aqj+i = a, 



where a is the initial and final spacing between two neighboring movers. We compute the 
so-called probability of reunion — the probability that the movers go back to their (almost 
coinciding) positions after time t. Of course, this is an exponentially vanishing probability. 
Notice that, since the A's also label highest weights of irreducible representations of U{N), this 
boundary condition is labeled by the Weyl vector for a suitable choice of the overall scale. Now 
a straightforward computation yields 

1 ^ 

PtA^o, Ao) = j—y^, n(l - e^'^yr-' . (5) 

We are going to relate this probability to the partition function of Chern-Simons theory. 

Recall that Chern-Simons theory is a topological quantum field theory whose action is built 
of a Chern-Simons term involving as gauge field a gauge connection associated to a group G on 
a three-manifold M [12] (see [15] for a recent review). The action is: 

SiA) = ^ f Tr(AAdA + -AAAAA) , (6) 

4-K Jm \ 3 / 

with k an integer number. Now if we choose units where a? = 1 and identify 

-t='' = kTN^ ^^^ 

with gs the string coupling, equation (5) is the partition function of Chern-Simons on S^ with 
gauge group U{N) [12]. Observables in Chern-Simons theory always come with a choice of 
framing, corresponding to a choice of trivialization of the tangent bundle in the gravitational 
Chern-Simons term [12]. In our case, notice that the framing is the matrix model framing 
[16, 17], which is related to the canonical framing as follows: 

Zcs{S') = ef'^\-^^^(^'-'^PtM>^o,Xo) , (8) 

where the label refers to the Weyl vector Aq = p, and as usual q = e^" . 

One way to understand why we get the partition function of Chern-Simons theory is to 
generalize the above to other compact groups [11, 18]. Using the method of images, one finds 
that the above probability generalizes to 

P*,r (A, n) = -^ e^^^ Yl <^)e^^''"''^'\ (9) 

where r is the rank of G and W the Weyl group. From here, using the Weyl denominator 
formula we can get amplitudes for more general boundary conditions: 



/^ ^ 1 p-P+IpI^ -tt (a, A) , , 

P^A\P) = (^;^^e-^^ n^-hi^ , (10) 



where a are the positive roots and p is the Weyl vector. This expression is the (unnormalized) 
expectation value of a Wilson loop around the unknot. The partition function is obtained by 
setting X = p. Normalizing (10) by the partition function gives the quantum dimension. 

Somewhat more generally, taking into account the framing and the central charge we can in 
fact write, dropping an overall sign, 

Pi,,(A,^)=e'^'i"-5(r5T)A^, (11) 



where 



I 2 > e{w)e ' 
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(X,w-fi) 



Txf, = Sx^e^(^^ 24 . (12) 

The central charge is c = kdiicng/{k + g), C(A) is the Casimir of the representation A, A+ is 
the set of positive roots, P is the weight lattice, and Q'^ is the coroot lattice. Obviously, S 
is the Brownian motion probability (with the external factors of T amputated), and T is the 
Boltzmann factor. It is now also clear that p{X, fi) itself corresponds to the (unnormalized) 
expectation value of the Hopf link invariant with representations A and /U. TST is in fact the 
operator that performs the modular transformation r -^ t/{t + 1) (see also [16]). By this 
particular surgery [12] one obtains S"^ out of two solid tori. 

To end this section, let us remark that there is a matrix model expression for the partition 
function of Chern-Simons on S"^ [19, 17]: 



..,.3,.£:i!pl!l/ni^eH.V..n(--^)^ . (13, 



1=1 i<j 

It is not hard to check that this is nothing but the extensivity property of probabilities: 



Pt+t',r{p,p) = [<iX\pt,riP^^)Pt'A^^P)^ (14) 

where the range of integration is the same as in the matrix model, but using symmetry can be 
restricted to the Weyl chamber. This can also be seen as a renormalization group property. It 
is also clear that the repulsive sinh, coming from the Vandermonde determinant, is related to 
the fact that the walls of the Weyl chamber are effectively repeling. Indeed, all the paths that 
end on these walls are suppressed from the expression for the probability, and hence only paths 
will contribute for which the particle stays away from the walls. 

We saw that the partition function comes out in the natural matrix model framing. It would 
be interesting to see if one can obtain more generic framings by rescalings and shifts of the 
boundary conditions. For the case of the partition function and the unknot this seems possible 
[8]. This results in real exponential factors, as expected from the analytic continuation of the 
phase factors. 

3 Brownian motion and QCD2 

It is now natural to ask what happens if we consider Brownian motion in the fundamental Weyl 
chamber of an affine Lie algebra. The affine Weyl group isW = Wt<T, where T is the group of 
translations in root space. In the case analyzed in the previous section, constraining the motion 
to the fundamental Weyl chamber was achieved by appropriately adding all images generated 
by the action of W. Now, in addition, we have to mod out by translations in the coroot lattice. 
We get the following density 

^^A>^'^^) = 7^ E E ^He--'-^^--^ (15) 

*^ ^ -yeiQ"^ wew 



where r is the rank and I = k + g. Standard manipulations yield 

for affine vectors A = A + /wQ) following common notation ujq = (0; 1; 0). In fact, the normalized 
density is the affine character chc(/i/t), in complete analogy with the finite case. With our 
boundary conditions, formula (f6) agrees with the general result found in [20]. 

Using the results in [21], it is not hard to see that the above is related to the partition 
function of 2d Yang-Mills theory on the cylinder. Two-dimensional Yang-Mills (see [22] for a 
review) is not a topological theory, but it is close to that in the sense that it has no local degrees 
of freedom. It is invariant under area-preserving diffeomorphisms. On the cylinder, the partition 
function is [23]: 

^2dYM(ff,5';t)= E XA(5^')XA(5')e^^^^")-5l''l\ (17) 

AeP++ 

where -P++ are the dominant weights^ . The gauge coupling Qym and the area of the cylinder A 
always appear through the combination t = ^ym^- 

From the cylinder one can easily obtain the partition function for other two-dimensional 
manifolds. A salient feature of the partition function on the cylinder is its extensivity. If we 
glue two cylinders of areas Ai and A2 along a common boundary, the partition function on the 
resulting cylinder of area A' = A1+A2 will have the same form, with new area A' . An analogous 
property also holds on the disk, and can be seen as an invariance under renormalization group 
transformations and indicates that the plaquette model for this theory is in fact exact. Notice 
that the area is linearly related to the Brownian motion time t, and so this is like the additivity 
property (14) that we encountered before. 

As is well-known, the partition function of 2d Yang-Mills is a solution of the heat equation 
on a group manifold. It is the kernel of the heat equation defined such that it satisfies the heat 
equation for both g and g\ and in the limit t ^ it tends to a Dirac delta function 5{g — g') 
with respect to the Haar measure. From it one can construct the unique solution of the heat 
equation with specified boundary condition at t = [24]. 

Let us now state the precise relation between Brownian motion and 2d Yang-Mills. We 
consider the group elements g = e^'^*'*'/ , g' = e27n^/( Using the manipulations in [21], we find 

Z2,YM{9,9;t) = ^^(2vrzA/0l),(-2vrz/./0 ^^'^^^'^^ ' ^'^^ 

where qt^r is given in (16). The normalization factors Dp come from the normalization of the 
characters that enter .^2dYM) whereas qt^i is always unnormalized, as we also saw in the previous 
section. In particular, the normalization is independent of t and thus the t-dependence is totally 
contained in qt^r- The above relation is our main result in this section. 

It is well-known that the partition function of 2d YM on the sphere does not have nice 
modular transformation properties. The reason is that it is proportional to the derivative of a 
modular form, rather than the modular form itself. For the case of the torus, see [25, 26]. The 
modular properties of the partition function on the cylinder are easy to work out from the above. 
It suffices to realize that qt^ri^, ^) is an affine character, which is obtained from a theta function 
by summing over all images. We hope to come back to this issue [8]. Related work appeared in 

^For convenience we included a constant factor of |p| . 



[27]. Let us here remark that there is a special group element which one can insert, so that the 
partition function on the sphere does have nice modular properties. This is the special group 
element defined in [28, 24] such that the character has values 0, ±1. Let us denote by \a) the 
associated state. We can then interpret .^2dYM(fl;0 ^^ ^ partition function on the sphere with 
insertion of a state \a) at t = 0. The partition function can be written in terms of the Dedekind 
r/- function by Macdonald's 77-function formula [29]: 

^2dYM(a;t) = e-^<i^-^r/(t)'i^-« . (19) 

We will end this section with some comments on one-dimensional discrete random walks and 
their relation to growing Young tableaux and Wilson's plaquette model^. Brownian motion is 
a limiting case of this [11]. We will focus on the random-turns model [11, 13], where N movers 
are allowed to move on a (one-dimensional) lattice, but at each tick of the clock only one mover 
walks, and it can take a step left or right with equal probability. It is known [31, 13] that the 
probability of reunion after 2n steps in this model is given by the expectation value of 2n powers 
of a unitary matrix in a unitary ensemble of NxN matrices. This probability is also proportional 
to the probability of the largest increasing subsequence of a random permutation of n objects 
of having length < N. By the Schensted correspondence between random permutations and 
Young diagrams, this is also the probability for the top row of a Young diagram with n boxes to 
have length < N. Let us call this probability P{ln < N). The Gross- Witten model is obtained 
by studying the Poissonized quantity, i.e. studying all possible (Poisson distributed) Young 
tableaux whose top row has length Lx < N: 



00 ^ 
P{Lx<N) = e-^V-A"P(;„<iV) 

n=0 



e-^ 



E E ^A^, (20) 



n=0 fihn, fJ.i<N 

where Lx is the length of the top row of a set of diagrams that are Poisson distributed, and 
fi\- n means that fi partitions n, and d^ is given by the hook formula of the diagram. But -j is 
the Plancherel measure, which is also the probability to pick a Young diagram of shape /j, among 
a random set. Thus, the above is the grand canonical ensemble distribution for tableaux with 
top rows of lengths Lx < N, and log A is the chemical potential for adding a box anywhere in 
the diagram so that it remains a valid U{N) tableau. If we identify A with the gauge coupling, 
we obtain Wilson's lattice version of QCD2 

/■ 1 °° 1 1 

Zgw = / dC/exp[-TVF(C/ + C/t)] = e^P{Lx <N) = Y, WJ^ — ^2n(^i = J, A, = i) , (21) 

where A = 1/(7^. The trace is taken in the fundamental representation of U{N). The last equality 
relates the partition function to the random walks probability. The plaquette model was solved 
at finite N in [32]. In [33] it was found that the model has a third-order phase transition^ at 
large N . Thus, lattice QCD2 can be reformulated as a growing (or shrinking) Young tableau. 

^There are various possible reformulations of the discrete random walks problem, that we will not consider 
here: as a Hammersley process, a growing PNG droplet, etc. See e.g. [30]. 

^It is not hard to see that the phase transition precisely comes from the restriction on the number of boxes in 
the top row to be < A''. 



The lattice model of QCD2 is celebrated for its phase transition. This occurs at A = g^N = 2 
and is closely related to a similar phase transition in the probabilities of the random distribution. 
Indeed, [34] has proved a depoissonization lemma that bounds the value of the probability of the 
random distribution from above and from below with its Poissonized value. The phase transition 
for the random Young diagram occurs when the number of boxes grows as n ~ 2\/N . 

Finally, let us comment that also in the case of the growing Young tableau there is a limiting 
shape after appropriate rescaling with N . This limiting shape is given by a continuous (but non- 
differentiable) function, the discontinuity being of course at the point of the phase transition. 
This is a two-dimensional analogue of the limiting shape of the 3d partitions found in [2] . 

4 Discussion and outlook 

In this letter we have pointed out three connections: between Brownian motion in the fun- 
damental Weyl chamber of a simply-laced, compact Lie group, and Chern-Simons theory on 
S^ for the corresponding group; between Brownian motion in the Weyl chamber of an affine 
Lie group and 2d Yang-Mills theory; and between the random-turns model of discrete random 
walks, a two-dimensional melting crystal, and lattice QCD2. The fact that the connections are 
quite general — in the case of Chern-Simons, we get a statistical mechanical realization of the 
modular matrices — and work for various gauge groups, might lead one to think that there 
may be more to the relation between random/diffusion walks and gauge theories than just a 
mathematical coincidence, and one might hope to find physical applications. Therefore these 
connections immediately raise several questions. First of all, how far does the correspondence 
go? In the case of Chern-Simons, we saw that the agreement can be understood from the repre- 
sentation of the modular S- and T-matrices as non-intersecting Brownian motion probabilities 
and Boltzmann factors, respectively. Taking into account the role of SL{2,Z) in the surgery 
approach of Chern-Simons theory [9], one may hope that for manifolds other than S'^ at least 
the simplest cases may have an interpretation in terms of Brownian motion quantities. This is 
certainly worth exploring further. Also, the basic cases of the partition function, the unknot 
and the Hopf link can be easily obtained from Brownian motion. It would be interesting to see 
if more general knots can also be obtained. 

In the affine case, where one finds the partition function of 2d Yang-Mills on the cylinder, 
it would be interesting to see if one can extend the connection to expectation values of Wilson 
lines, or whether there is a Brownian motion interpretation of the partition function on the 
three-punctured sphere. 

A particularly interesting point would be to see if Brownian motion can give us further 
connections between these low-dimensional theories. For example, the fact that Brownian motion 
is a limit of a discrete random walk is very suggestive of a connection between two-dimensional 
and three-dimensional theories. Also, it would be interesting to see whether the lock-step model 
[11], which we have not considered in this paper, also has a reformulation in terms of a low- 
dimensional gauge theory. Notice also that the way 2d Yang-Mills arises from Brownian motion 
is by effectively compactifying the dual Cartan subalgebra to a torus. The resulting expression 
is an affine character, which makes the modular properties of the partition function completely 
explicit. A connection between 2d Yang-Mills on the torus and topological strings has recently 
been pointed out in [25]. 

We saw that the matrix model of Chern-Simons theory on S^ naturally arises in the compo- 
sition law of probabilities. Indeed, since we are dealing with continuous paths it seems that the 
matrix model formulation of Chern-Simons theory is the most natural one for Brownian motion. 



Yet from the point of view of the WZW model one naturally gets sums over representations 
rather than integrals. At the level of the intermediate states, there is a precise way in which 
both approaches are equivalent [8]. In Chern-Simons theory the representations are integrable, 
and from the Brownian motion point of view these correspond to special points on the line. We 
can deal more generally with arbitrary points by using characters. It would also be interesting 
to explore the connection with the fermionic representation. 

On the more mathematical side, the underlying principle allowing these connections seems 
to be the fact that all these models in one way or another satisfy the heat equation. It has been 
known for a long time that certain quantities in the WZW model satisfy the heat equation [35]. 
Also, the heat equation is closely related to modular invariance. We hope to report more on 
this in the future [8]. 

Perhaps one of the most interesting questions is whether Brownian motion can be used as 
a tool in string theory, in the spirit of [2, 4], for example. Chern-Simons theory is the effective 
gauge theory describing the topological A-model [9], and so a reformulation in terms of Brownian 
motion might be useful for string theory itself. Notice furthermore that the natural string theory 
coupling is related to the Brownian motion parameter t (which is actually the product of the time 
parameter and the diffusion coefficient) as —j = Qs, whereas the relation to the Chern-Simons 
coupling involves analytic continuation*^. This suggests that the interpretation of topological 
strings in terms of a statistical mechanical system may in some ways be more natural than as a 
gauge theory. In particular, it would be extremely interesting to understand whether the heat 
equation plays a role in the topological A-model. One would also like to see if discrete random 
walks — of which Brownian motion is a limit — are related to topological strings. Notice that 
[36] have used random partitions — which, as pointed out, are equivalent to the random walks 
model and, after Poissonization, to the plaquette model of QCD2 — to compute the prepotential 
oiN = 2 SYM theory. 

Another interesting question is whether 2d Yang- Mills and lattice QCD2, and their respective 
third-order phase transitions, have string theory interpretations. In [37] it was shown that the 
QCD2 plaquette model can be used to obtain the SU{2) J\f = 2 Seiberg-Witten solution by 
taking a double scaling limit. In particular, a local Calabi-Yau geometry that engineers this 
curve was found. In this case, it was argued that the phase transition plays no role. In [38], 2d 
Yang- Mills on a Riemann surface E was obtained by wrapping D6 branes on S*^ x S. 

We saw that the Gross- Witten phase transition does play a role in the context of the shrinking 
two-dimensional Young tableau. It was related to the non-differentiability of the limiting shape. 
It would be interesting to see whether such phase transitions are present and play any role for 
the topological vertex. 

Let us also mention that Cardy [39] has recently found a remarkable connection between a 
system of N non-intersecting Brownian motions (described through the celebrated SLE process) 
and boundary-bulk conformal field theory models and integrable models of Sutherland type. 

Finally, recall that [40] pointed out connections between vertex models and Chern-Simons 
theory. The random walks that we have mentioned in this paper are special cases of vertex 
models. However, for us the connection with Chern-Simons theory appears in the continuous 
limit of Brownian motion rather than in the discrete case. 



This analytic continuation is a subtle issue that we hope to get back to in the future. 
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